Asymmetric Rotational Stroke in Mouse Node Cilia during Left— Right Determination 
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Clockwise rotational movement of isolated single cilia in mice embryo was investigated in vivo. 
The movement generates leftward fluid flow in the node cavity and plays an important role in left- 
right determination. The leftward unidirectional flow results from tilting of the rotational axis of the 
cilium to the posterior side. Because of the no-slip boundary condition at the cell surface, the upper 
stroke away from the boundary generates leftward flow, and the lower stroke close to the boundary 
generates slower rightward flow. By combining computational fluid dynamics with experimental 
observations, we demonstrate that the leftward stroke can be more effective than expected for cases 
in which cilia tilting alone is considered with the no-slip condition under constant driving force. Our 
results suggest that the driving force is asymmetric and that it is determined by the tilting angle 
and open angle of the rotating cilia. Specifically, it is maximized in the leftward stroke when the 
cilia move far from the boundary. 
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Cilia or flagella are widely found in biological systems 
from bacteria to mammals [1]. Their movements relate 
directly or indirectly to their own biological functions, 
e.g., protozoa swimming, coordinated beating of cilia in 
epithelium, and rotational movement of mono-cilia in 
mouse nodes. This study specifically examined the third, 
mouse node cilia, which are involved in left-right (LR) 
determination 0, 

A small cavity, called a node, is formed in a mouse 
embryo at an early developmental stage (7.75 days after 
fertilization). The node cavity, which is filled with extra 
embryonic fluid, is constructed with a few hundred cells, 
each of which protrudes a single motile cilium (Fig. [Ha)). 
The cilium consists of nine doublet microtubules without 
a central pair, called 9-f or a primary cilium, differently 
from the other ordinary motile cilia consisting of nine 
doublets with central pairs of microtubules (9-1-2). Cilia 
generate a rotational or beating movement (mainly mov- 
ing back and forth), respectively, in 9+0 or 9+2 cilia, by 
sliding motion of dynein motor proteins along the micro- 
tubules. It is particularly interesting that the rotational 
movement is almost always clockwise when viewed from 
the ventral side, at least in mouse, rabbit, and medaka 
In this paper, we show that the driving force of the 
rotational movement can be rotationally asymmetric in 
contrast to the fact the 9+0 has structural rotational 
symmetry. 

The clockwise rotational movement by the 9+0 cilia 
serves an important role in LR symmetry breaking. The 
movement generates leftward laminar flow of up to 1- 
10 /im/s because the rotational axis of the cilium tilts 
to posterior side. Because of a no-slip boundary (NSB) 
condition on the cell surface, upper leftward strokes (LS), 
which are less affected by NSB, generate faster leftward 
flow, and a lower rightward stroke (RS), which is directly 
affected by the boundary, generates a slower rightward 
flow (Fig. mb)). Results reported herein show that LS is 
more effective than expected before, when only the NSB 



condition was considered under a constant driving force 
of the rotational movement 2']. The result was derived 
by analyzing rotational movement of isolated single cilia 
using mutant mice, then combined with computational 
fluid dynamics. 




FIG. 1. Node cilia. L, R, and P respectively denote left, right, 
and posterior directions in a node: NC, node cell; NF, node 
fiow; LS, leftward stroke; RS, rightward stroke; BB, basal 
body; IJ, intercellular junction. 



Rotational movement of a single isolated cilium ob- 
tained from a ciliary mutant Dpcd Q was recorded using 
high-speed imaging for 2 s (100 frames/s; HAS-SOOM; Di- 
tect Co. Ltd.) 6]. The cilia tip portions were tracked 
manually to obtain 2-dim projection on the x-y plane 
(Fig. mja), Figs. EJa-c)). The cilia setup and movement 
were characterized using the following parameters; the 
tilt angle of the rotation axis a, the open angle j3 of the 
rotation orbit along a conical face, the tilting direction 77, 
cilium length Z, and rotation phase 9 (Fig. [2]). The phase 
9 is measured from the position closest to the cell surface. 
The projected trajectory to the observation plane {x-y) 
is elliptic as a result of tilting of the cilium. Therefore, 
9 was estimated as that after applying ellipse-correction 
[41 using multivariate nonlinear regression. Radii of the 
major and minor axes a and 6 were obtained from the 
fitted ellipse. The tilting direction 77 was calculated from 
the minor axis of h combining with the root position of 
cilia R. The other parameters were calculated with the 
following relations derived from the construction shown 
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x-y plane 



FIG. 2. (color online) (a) Setup of cilia, (b) top view (x-y 
plane projection), (c) side view (z-ry projection), R denotes 
the cilia root position. 



in Fig. mb, c): a = cos ^(fe/a), I = a^l + ^^^b^, 

/3 = sin^^(a/Z), where s denotes the distance between 
R and the center of the eUipse. Estimation of a includes 
some uncertainty (roughly assumed to be ±10° (degree) 
= 0.17 (rad)) because of experimental setup problems: it 
is difficult to put the node samples perfectly fiat against 
the observation plane. 

To infer the driving force from the observed data, we 
used computer-aided fluid dynamics, where a cilium is 
modeled as a rigid cylinder with length of I and diameter 
of r rotating around a fixed orbit. The driving force is 
balanced with the sum of the hydrodynamic drag force 
and the elastic restoring force. The magnitude of the 
elastic force must be almost zero or constant because cur- 
vature of the cilium is unchanged during the rotation as 
observed in the experiment 0. Additionally the normal 
component of the drag force is much smaller than the tan- 
gential component. Consequently, the driving force can 
be approximated as the one balanced by the tangential 
component of the drag force. The equation of rotational 
movement in the Stokes flow regime is written as 



m 



(1) 



where F{0) is the driving force required for the rota- 
tional motion, and J{9) is the drag coefficient which is 
calculated numerically using boundary element method 
(BEM) as follows 0- 

When a cilium moves in a fluid bounded by an infinite 
flat wall (approximation for a surface consisted of an ar- 
ray of cells) , the flow field can be given as the following 
integral form 

u(v) = --^ / G(v - w)q(w)dA, (2) 

the vector v is a point where the velocity u is calculated, 
and w is a point where the traction force q is exerted. 



/i represents the viscosity, and A denotes the cilium sur- 
face. The function G is the half-space Green function 
[Tol |. which rigorously satisfies the wall boundary condi- 
tion [llj. No background flow of the surrounding fluid 
is assumed. The NSB condition on the wall as well as 
the surface of cilia was assumed. The fluid has the same 
velocity to the rotating cilium on its surface. Using Eq. 
([2]), one can calculate the traction force distribution on 
the cilium surface from the velocity distribution. Equa- 
tion ([2]) is solved numerically using BEM |12| , where each 
cilium surface is discretized with 594 triangular elements 
(Fig. [3l^a)). The integration in Eq. ([2]) is performed on 
a triangle element using 28-point Gaussian polynomials. 
The singularity in the integration is solved analytically 
The flnal linear algebraic equation is solved by LU 
decomposition, including pivoting. Because of the linear- 
ity of the Stokes flow, the drag force is proportional to the 
rotational velocity 9 [l^. Consequently, the strength of 
the drag force Fd becomes equivalent to J {9) when ^ = 1, 
where Fd is deflned as Fd — — St • qdA. Therein, St 
is the unit tangential vector along the circular trajectory 
of the top of the cilium. Once Eq. ([2]) is solved for q, Fd 
can be calculated, so does J {9). In the numerical calcula- 
tions, the parameters characterizing each cilium obtained 
from the experimental observations described above were 
used. 

Figure Efb) presents examples of numerical results of 
J (9) . The drag force is maximized when the cilium passes 
closest to the cell surface (9 = 0). The mean value and 
the amplitude of the drag force increase when a, /3, or 
I is enlarged. The mean values of J{9) range 0.05-0.10 
(pNs) as shown in Figs. \M,c) andlHt^g-i). These values are 
slightly larger than those calculated from the relation on 

a drag coefficient i^fjy-^'^Q 447 [l^ ■ which is conventionally 
used for torque estimation of rotating molecular motors 
In this relation, the NSB condition is considered, 
but a cilium is set parallel to the boundary. Therefore, 
tilting of the cilium is not considered. 

Now the numerical results of drag coefficients J {9), and 
experimental observations of 9 and 9 are combined by Eq. 
(d]). Using formulations F{9) = Fo(l - FrC0s{9 - F,)) 
and J{9 — Js) ( replacement of J{9) with phase shift Js 
), three models are proposed as shown below. 

(1) constant F {Fr = 0, = 0), 

(2) constant F and shifted J (F,. = 0, Js ^ 0) 

(3) variable F depending on 9 {0 < Fr < 1, Js = 0) 

Here, phase shift Js was introduced to assess the pos- 
sibility of estimation error of 77, which might affect the 
definition for the position of 6* = 0. Parameters Fq, Fr, 
and Fs represent the mean driving force, force skcwness 
ratio, and phase shift of force maximum from 9 ~ tt (the 
phase a cilium pass a position furthest from the cell sur- 
face), thereby Fq > 0, < Fr < I, and -tt < < tt. 
Figure m^b, c) shows one application. The function J{9) 
was calculated using BEM based on the ex per imental ob- 
servation of a, /3, / and r = 0.4 {fim) [l3| |ll |. 
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FIG. 3. (color online) (a) Setup for boundary element 
method. Discretized cilium as a rigid cylinder (red), (b) Drag 
coefficient J (6) normalized into dimensionless form by fj,ujl^, 
where jj. is the viscosity, and where uj is the given rotational 
speed. The thick solid line shows the result obtained when the 
parameters was set as Z = 1, a = n/6, /? — it/4:, r — 0.1. Thin, 
dotted, and dashed lines respectively represent results when 
one parameter changed to I — 1.5, a = 2n/9, and /3 = 2ir/9. 
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FIG. 4. (color online) Example of estimation of F{9). (a) 
Projected trajectory, where cilia parameters were estimated 
as Q = 0.37, 13 = 0.61, I = 3.2 (/im). (b) Phase speed fitted 
by Eq. ((l| using the three models. Green dotted and dashed, 
blue dashed, and red lines respectively represent fitting results 
by model (1-3). Data numbers of phase-speed n = 44. (c) 
F{9) calculated using model (3). Black solid lines denote 
J (6), red dashed lines denote F{6). Red short dashed lines 
represent a 95%- confidence interval associated with F{6), 
which is estimated from the dataset of phase and phase speed 
(plots of Fig. IHb)) using the bootstrap method ila |. 



To consult the degree of fitness statistically, model se- 
lection analysis was applied using the Akaike Information 
Criterion (AIC) The AIC value is a measure of the 
relative goodness of fit of a model, defined as —2 \nL+2k, 
where L is the maximum likelihood function and k is the 
number of parameters. Therefore, a smaller value of AIC 
refiects better fit of the model. Table U presents a sum- 
mary of the result for the case of Fig. The estimation 
on a includes some uncertainty. Therefore, three vari- 
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FIG. 5. (color online) Calculation examples of F{6). (a-c) 
Projected trajectory. Cilia parameters were estimated respec- 
tively as Q=0.30, 0.10, 0.68, ^=0.58, 1.24, 0.86, I =3.4, 2.4, 
2.8 {fJ,in). (d-f) Phase-speed fitted by the three models, (g- 
i) ^'(6') calculated using model (3). Other notations are the 
same as those of Fig. |4] 



ations (estimated a ± 0.17) were tested. Consequently, 
model (3) when a = 0.37 (No. 1) shows the best fit. Fig- 
ure [5] shows other examples, and model (3) best explains 
the experimental data in any case |20j |. 



TABLE I. Model selection: AIC value for each model (1-3). 
The other parameters except a are the same as those in Fig. 
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No. 


Q:(rad) 


model (1) 


model (2) 


model (3) 


1 


0.37 


300.2 


300.4 


275.8 


2 


0.54 


293.7 


292.1 


277.2 


3 


0.72 


287.7 


280.7 


278.2 



Figure ini summarizes estimated parameters Fr, Fs, and 
Fq for F{0) depending ona+P [2l|, of which the comple- 
mentary angles correspond to the distance between a cil- 
ium and cell surface when the cilium passes closest to the 
cell surface. The force skewness ratio Fr is 0.1-0. 5s (Fig. 
[6fa)), suggesting the existence of asymmetric stroke that 
cannot be explained only by the effect of tilting cilia with 
the NSB condition at the cell surface if a rigid straight 
cylinder for cilium can be assumed. Straightforwardly 
interpreting this result, the driving force is not constant 
and maximized at a phase 9 = n + Fg (Fig. IHb)). In 
other words, it is maximized in LS and minimized in RS 
when a + l3 is small (up to ~ 37r/8) as shown in Fig. [lUd), 
whereas it is maximized at the first half of RS when a + jS 
close to tt/2 as shown in Fig. HJe). 

The mean driving force Fq was around 2 (pN) except 
when a+P is close to 7r/2. The order of magnitude is 
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FIG. 6. (color online) (a-c) Relation between a + /3 and pa- 
rameters of driving force F(d) — Fo{l — FrCos{6 — Fs)) es- 
timated using model (3) from 11 trials obtained from four 
isolated cilia. Vertical and horizontal error bars represent 
a 95%- confidence interval associated, respectively, with the 
parameters of F and with a + fi, which is estimated using 
the bootstrap method fl3], respectively, with the dataset of 
phase and phase speed, and ellipse fitting. Solid lines repre- 
sent fitting by appropriate curves: ax'' (a, b) and a/{b — x) + c 
(c). The parameters were estimated using nonlinear regres- 
sion: a = 0.19 ± 0.05, b = 1.03 ± 0.73 (a); a = 0.57 ± 0.24, 
b = 3.68 ± 1.14 (b); a = 0.096 ± 0.048 b = 1.57 ± 0.01, 
c — 2.01 ± 0.23 (c). (d,e) Schematic diagram for asymmetric 
driving force when q + /? is small (d) and large (e). 



equivalent to the force generated using a single dynein 
head sliding along a microtubule in 9-1-2 cilia or flagella 
(estimated as 2-5pN in the literature) [11,[2^, which sug- 
gests that a single dynein participates in the generation 
of rotational movement at a moment. 

The value Fq grows rapidly when a+P approaches 7r/2. 
The cilium skimmed over the cell surface in RS, then 
F{9) maximized at around 3tt/2 < 9 < 2tt (Fig. El^e)). 



This can be understood from the observation that force 
generation of 9-1-2 cilia depends on a load jlH- This 
fact, however, does not explain how cilia accelerate in 
LS more than in the case expected by tilting cilia with 
NSB condition when a+P is much smaller than 7r/2 (Fig. 
ini[d)). The following four ideas are conjectured as an- 
swers. First, rotational motion is regarded as generated 
by a dynein, protruding from one of nine doublets of mi- 
crotubules attaching the next doublets sequentially in a 
circumferential direction in nine doubles [23j] . The phase 
speed is greater in LS even if only tilting cilia with the 
NSB condition are considered, which can raise the prob- 
ability of a dynein attaching to the next doublet. This 
attachment can accelerate the rotation speed in LS and 
vice versa in RS. Second, tensions acting on certain cy- 
toskeleton connected to a basal body can be asymmetric 
[2^ 1 because of shifting of the basal-body position in each 
cell (Fig. [TJb)) [2^. Third, axonemes themselves can 
be asymmetric. Such a structure is found in 9-f 2 cilia 
(26j . Finally, the arcuate shape of cilia can also affect 
the asymmetric stroke p7| . 

In conclusion, our results suggest the existence of 
asymmetric strokes that cannot be explained adequately 
by the effect of cilia tilting with NSB condition at the 
cell surface. This would result from the asymmetric force 
generation in LR strokes if a rigid straight cylinder for 
cilium can be assumed. In an early stage of develop- 
ment, the number of cilia, length, rotational speed, and 
tilting angles are small [1, [25|. However, it was revealed 
recently that sufficient flow speed for LR determination 
can be achieved by only two cilia @. Under such an 
unreliable situation by a few cilia, the asymmetric stroke 
can support the establishment of LR determination more 
effectively. 
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